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Abstract 



The Characteristic Angle (CA) method [Lei Zhou and Ruibao Tao, J. Phys. A 
27 5599] developed previously for the easy-plane spin-1 magnetic systems has been 
successfully extended to the spin-| case. A compact form of the CA spin-| operator 
transformation is given, then the ground state energy, the magnon dispersion relation 
and the spontaneous magnetization are discussed for an easy-plane spin-| ferromagnet 
by using the CA method. Comparisons with the old theoretical methods are made in 
the end. 
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1 Introduction 



"Easy-plane" single-ion anisotropy (D(Sf) 2 ,D > 0) is very popular in magnetic sys- 
tems with spin greater than one half, and such systems have been attracting much at- 
tention for years [1-11]. On theoretical side, the spin-wave excitations in such systems 
are very difficult to handle caused by the off-diagonal terms introduced by the "easy- 
plane" single-ion anisotropy term. Actually, many straightforward methods which 
are good to deal with the "easy-axis" single-ion anisotropy (D(Sf) 2 ,D < 0) are usu- 
ally found to be invalid for the "easy-plane" magnetic systems no matter how weak 
the anisotropy parameter D will be [4,11-13]. For example, applying the Holstein- 
Primakoff (H-P) transformation [12] naively to study the magnon excitations in the 
easy-plane magnetic systems must lead to an imaginary value for the energy of the 
"k=0" mode [4,11], and introducing a vector rotation of the spins in the XZ plane 
to optimize the spontaneous magnetized direction is also helpless [13] - the excitation 
energy in a harmonic approximation under such an approach will either be negative 
or be imaginary. In fact, caused by the off-diagonal terms (D((Sf) 2 + (S^) 2 )) in the 
Hamiltonian introduced by the single-ion anisotropic terms, the proper eigenstate of 
the Hamiltonian must be a mixture of the single-site spin states \n) with \n + 2) and 
\n — 2) [3,4]. Such a spin-states mixing effect is very significant in the present "easy- 
plane" the result, an ordinary method which fails to consider such an effect 
must be invalid for such systems. The matching of matrix elements (MME) method 
is one possible way to take this effect into account perturbatively [3,4], and many ap- 
plications of this method have been implemented in a first order approximation of the 
anisotropy term [3-7]. Some numerical methods had also been developed for the easy- 
plane spin-1 ferromagnetic systems where the authors had considered the spin-states 
mixing effect [9-10]. Recently, a new method - the characteristic angle (CA) method 
was developed successfully for the easy-plane spin-1 magnetic systems [1], in which the 
spin operators are transformed to a new set of quasi-spin operators which have taken 
account of the spin-states mixing effect by a variation parameter - the characteristic 
angle (CA). The ground state properties of an easy-plane spin-1 ferromagnet, such as 
the ground state energy, the spontaneous magnetization and the magnon dispersion 
relation, are more reasonable in the CA approach than those in the MME method [1]. 
The induced magnetization of the easy-plane spin-1 ferromagnet has also been studied 
based on such an approach [11]. 

In this paper, we would like to extend the CA method to the easy-plane spin-| 
magnetic systems. Although the main idea has been proposed in Ref. 1, an extension 
of this approach to larger spin case is still nontrivial because the extension makes it 
possible to study the magnetic systems with single-ion anisotropy in larger spin case. 
Actually, the extension is not very easy because the number of the variation parameters 
will increase and the expression of the CA transformation will become more complicated 
than that in spin-1 case. This paper is organized as follows. After diagonalizing the 
single-site Hamiltonian part, we derive the CA spin operator transformation for S — | 
case successfully in Section 2, then we present the main physical results for an easy- 
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plane spin-| ferromagnet using the CA method in Section 3, and the comparisons with 
the MME method and other methods are made thereafter. Finally, conclusions are 
outlined in Section 4. 



2 Characteristic Angle transformation for S = | 

The Hamiltonian of an easy-plane spin-| ferromagnet can be given as 



(1) 



where (i, j) means summation restricted on the nearest-neighbor pairs. The spin vec- 
tors are forced into the YZ plane by the anisotropy term (D > 0), so that the sponta- 
neous magnetized direction must be in the YZ plane. Without loosing any generality, 
the Z direction is chosen to be the spontaneous magnetized direction and an external 
magnetic field h will be applied along the Z direction which will be set to zero in the 
end. 

Following Ref. 1, let us start from diagonalizing the single-site part of the Hamil- 
tonian. 

The matrix form of the single-site Hamiltonian D(S X ) 2 — hS z in the S z representa- 
tion is found to be 
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where the matrix elements are defined by A mn = (m\D(S x ) 2 — hS z \n) in which 
the bases \m) are the complete set of eigenstates of operator S z with eigenvalues 

V2' 2' 2' 2'' 

A new representation can be defined as follows, in which the complete set of bases 
\rh) are related to the old ones by the following orthogonal transformation: 
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Then, after the transformation, the matrix form of the operator D(S X ) 2 — hS z in 
the new representation \fh) will be 
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A = V T AV. 
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It is easy to check that only when 
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the matrix A can be exactly diagonalized with the following eigenvalues 

lh - VD 2 + h 2 + Dh 
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Since the orthogonal transformation V can exactly diagonalize the single-site part 
of Hamiltonian, then if we apply this transformation to the total Hamiltonian (1), we 
can expect to get a reasonable representation of the easy-plane spin-| ferromagnet after 
selecting appropriate values of the parameters 9\, 9 2 by the variation method. 

However, for the sake of simplification in this paper, we will only study the easy- 
plane spin-| ferromagnet in zero external magnetic field here (in another word, h — 0), 
then the two variation parameters 9i and 9 2 should be the same. In fact, since the 
Z direction is not specified in Hamiltonian (1), the state that all spins are at ||) and 
the state that all spins are at | — §) will give the same physical results. So, from the 
definition of the transformation V (Eq. flD), we may easily obtain 9\ = 9 2 = 9. Of 
cause, if an external magnetic field h is applied, the two parameters may not be the 
same and things are somewhat complicated then. 

In order to investigate how the transformation acts on the exchange interactions, 
we should first study how the transformation acts on the spin operators. The matrix 
form of the spin operator S + in the old S z representation is 
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After the orthogonal transformation, the matrix form of operator S + in the new 
representation is found to be 
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Now we may understand this transformation in another way: suppose the complete 
bases \m) are maintained unchanged, instead the spin operators are transformed to 
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new ones by an unitary transformation. Noticing some matrix identities 
operators such as 
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we can easily get the following spin-operator transformation: 
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Eqs. ( fL3|) -(|l5D are just the Characteristic Angle (CA) spin operator transformation 
for the spin-| case, and 9 is the Characteristic Angle (CA) which is actually a variation 
parameter determined later by minimizing the ground state energy. The transformed 
spin operators, so called CA spin operators (S ± ,S Z ), obey all the spin-| operator's 
commutation rules. The proof of the above CA transformation is straightforward. One 
may only write out the matrix form of the CA spin operator S + (Eq. (|T3"D ) in the 
ordinary \m) representation, then can find out it is the same as the matrix on the right 
side of Eq. (|i~0|) . This means the CA transformation for the spin operators is actually 
equivalent to the orthogonal transformation for the spin space. Furthermore, since an 
orthogonal transformation does not affect the commutation relations for matrices, and 
since the matrix form of the CA spin operator is related to the matrix form of the 
original spin operator by an orthogonal transformation, the CA spin operator (S^, S z ) 
must obey all the spin-| operators' commutation rules. 

By the way, although it seems possible to obtain the CA transformations for arbi- 
trary spin case following the method described in this section, however, a closed form of 
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the CA transformation is not easy to derive for larger spin case because more variation 
parameters will be needed then. An extension of this method to the general case is 
very difficult. 



3 Results and Comparisons 



In order to study the ground state properties and the low-lying excitations for such a 
system, the H-P transformation is introduced for the spin-| operators: 
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Applying the above H-P transformation, we can obtain the Bose expansions for CA 
spin- 1 operators. 

However, there exist another straightforward method to determine the first several 
order terms of the Bose expansions for the CA spin operators. We will consider S + 
as an example. After transforming the Hilbert space of spin operator to that of the 
Bose operator, the vacuum state of Bose operator |0)b is defined to be the spin state 
1 1)5, and the excited Bose states |2)b, |3)b are the corresponding spin states: 

\\) s,\ — \) Si\ — \) Si respectively (The subscript B and S are used to denote the Bose state 
and the spin state, respectively). Other Bose states such as \n) B , n > 3 are unphysical 
states. In physical space, the matrix forms of some lower order Bose expansion can be 
written out readily: 
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etc. Then from the matrix form of the operator S + (Eq. (]TD|) ) , it is very easy to write 
out the first order terms of CA spin operator S + : 



cos 29 a — sin 29 a + + 
The Bose expansion of S z can be similarly derived as: 
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Other terms can be obtained easily. 

Applying the CA transformation to Hamiltonian (1), we get 

H = -JY,^^i + DY,&) 2 - (23) 

(hi) 1 

The transformed Hamiltonian H must have the same eigenvalues with the original one 
since the transformation is an orthogonal one. 

Then we can apply the H-P transformation to the Hamiltonian Eq. (p3|) keeping a 
harmonic approximation. We get 

H = U + H 2 + --- (24) 
U = [- JZ^- - 2 sin 2 0) 2 + - D cos 2 6 - ^j-Dsm26)]N, (25) 

H 2 = JE{( 3 - 4sin2 ^)[ a ^-^ sin2 ^(K + ) 2 + K) 2 )] 

- [(3 cos 2 26 + sin 2 26)af aj + V3 sin 29 cos 26(afaj + a^)]} 
+ D Ei^Oy cos 261 - i sin 26) [(a+) 2 + (a,) 2 ] 

i 

+ (V3sin20 + cos20) a+a^}. (26) 

Taking a Fourier transformation to the momentum k space, we get the following 
Hamiltonian: 

H = U + J2 A kata k + J2 B k(atat k + a k a_ k )^ (27) 

k k 

A k = JZ(3-4sin 2 #) - JZ(3 cos 2 20 + sin 2 2#) 7fc 

+ v / 3£'sin2# + J Dcos2#, (28) 

B k = — [-JZ(3 -4sin 2 #)sin26) + — Dcos26- -Dsin26] 

+JZV3sm26cos26-f k , (29) 

where 

lk = ^^exp(ik-r), (30) 
z 8 

in which the 5 summation runs over the Z nearest-neighbor sites of a given site. 

The approximate Hamiltonian Eq. fl2~7|) can be diagonalized by a Bogolyubov trans- 
formation, 

H = E Q + Y J ^ata k + --- (31) 

k 

Eo = U,-W^ + \Y.^l-^Bl (32) 

k k 
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The spontaneous magnetization for the easy-plane spin-| ferromagnet in the harmonic 
approximation can be derived as: 

i ]-A k + y/2sm26B k 
M = (0\S Z \0) = 2-2sm 2 6 - —Y 2 . (34) 

As we know, the value of variation parameter 9 should be determined by minimizing 
th ground state energy E . An analytical expression of 9 with respect to the anisotropy 
parameter seems very difficult to obtain, however, an expansion in d = D / JZ can be 
determined order by order. As a first step, suppose d is very small so that 9 is also 
small, then by minimizing the ground state energy keeping only the first order terms 
in d, we obtain the following equation: 

—E~—U = —2JZ sin 49 — (2JZ + D) sin 29 + y/3D cos 29 
d9 d9 

~ -2.1 Z ■ 49 - 2.1 Z ■ 29 + V3D = 0, (35) 

then, 

9 ~ —d. (36) 
12 v ; 



Applying this solution into Eqs. (|2l|)-(|22T), it is easy to obtain 

S+ = s /3a-—da + --- (37) 
6 

S~ = V3a + -—da--- (38) 
6 

S z = ^- a +a + ^d{a +2 + a 2 ) ■ ■ ■ , (39) 

then from Eqs. fl3"2"|)-(|3"3]), we get the ground state energy and the magnon dispersion 
relation as following: 

Eo = ~JZN(^ - id) - ~ JZJ2 M 1 ~ Th) + d ? ~ (« 2 , (40) 

k 



e k = JZ^[?,{l- lk )+df-{d lk f. (41) 

It is easy to check that e k — > when k — ► 0. 

Comparing the transformation Eqs. (|37j)-(|3"9|) and the results Eqs. (^0|)-(^T]) to 
those of the MME method which are obtained in a first order approximation of d [4,6], 
we are interested to find that they are exactly the same. 

In the case of large anisotropy d, the first order approximation is not enough, numer- 
ical calculation must be required. We have studied a simple-cubic lattice in the present 
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paper, the characteristic angle 6, the ground state energy Eq and the spontaneous mag- 
netization M have been shown in figures 1-3 as functions of the anisotropy parameter 
D/JZ. The most interesting thing is that the easy-plane spin-| ferromagnet does not 
exhibit a phase transition from the ferromagnetic phase to the non-ferromagnetic one 
as D/JZ is enlarged (figure 3). This is quite different from the spin-1 systems [1,9]. 
Actually, when D ^> JZ, the single-ion term is dominant. In this case, one may find 
from Eqs. ©-(0) that 9 = tt/6. According to Eq. (|22), we get: when D ^> JZ, 
M — ► 1. However, in spin-1 systems, one may find that when D > D c (< 4JZ), 
M — ► [1,9]. 



In the remaining part of this section, we will compare the CA method with the 
MME method and other old methods in details. The total Hamiltonian of an easy- 
plane spin- 1 ferromagnet can be divided into two parts: 

H = -jY J {SlS z ] +StSj) + DY J {S*) 2 = H Q + H mt + const. (42) 
H = -2JZSY,S? + DY,m 2 , (43) 

i i 

H mt = -Jj2[StSy + (S*-S)(S*-S)}. (44) 

The main question in such systems is to find a best "starting point" to perform the 
spin waves calculation. It is just because the "starting point" are chosen unreasonably 
that the methods mentioned in the introduction are invalid for the easy-plane magnetic 
systems [4,11-13]. For example, the "starting point" in H-P transformation method 
is the ordinary ferromagnetic state that all spins are aligned parallel to each other, 
however, this "starting point" does not consider the spin-states mixing effect at all so 
that it must be unreasonable. 

According to the MME method [4], the "starting point" is the ferromagnetic state in 
the diagonalized representation of Hamiltonian H (Eq. (f£3|)). Treating the anisotropy 
term D(S^) 2 as a perturbation, the Hamiltonian H can be diagonalized in terms of 
d = D/JZ order by order. Based on this, the spin-Bose transformation was applied 
to the total Hamiltonian to study the magnetic properties of the system. In such an 
approach, it should be noted that the Hi nt part of the Hamiltonian has been neglected 
when determining the "starting point" , so this approximation may be valid only when 
the anisotropy is very small. 

However, we believe that the "starting point" in the CA approach will be more 
reasonable because the H int has also been taken into account when determining the 
"starting point" . As we can see in this paper, the CA transformation is applied to the 
total Hamiltonian and 9 is determined by minimizing the total ground state energy 
Eq (Eq. (|3"2"D). Actually, H int has been considered through its contribution to Eq. 
Since H int will contribute the second order terms in d to the ground state energy E , 
it is very easy to understand that the CA method gives the same physical results as 
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the MME method in the first order approximation of d which has been demonstrated 
in the present work. However, for large d case when higher order terms in d should 
be considered, the CA method may give more accurate results than the MME method 
because the CA method has considered more contribution than the MME method when 
choosing the "starting point" . 

Of cause, the CA method can only deal with the spin-1 and the spin-| cases at 
present time, however, the MME method can deal with arbitrary spin case perturba- 
tively. 

4 Conclusions 

In this paper, we have successfully extended the characteristic angle method to the easy- 
plane spin- 1 magnetic systems. Based on the diagonalization of the single-site part of 
the Hamiltonian, a compact form of the CA transformation for the spin-| operators 
has been derived, in which the spin-states mixing effect is considered automatically by 
the variation parameters. After applying the CA transformation and the H-P spin- 
Bose transformation, a model Hamiltonian which has the same eigenvalues with the 
original Hamiltonian is presented by the Bose operators. Bogolyubov transformation is 
applied to diagonalize the Hamiltonian in harmonic approximation, and the variation 
parameter CA is determined by minimizing the ground state energy. Analytical results 
of the magnetic properties are obtained in a first order approximation of d = D/JZ 
which are found to recover the MME method. Numerical calculations are carried out in 
large d case for the easy-plane spin-| ferromagnet, and detailed comparisons with the 
MME method and other methods are made in the end. We find that the CA method 
has considered more contributions than the MME method and other methods for the 
spin-1 and spin-| easy-plane magnetic systems. 
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Captions: 



Figure 1: Value of the characteristic angle as the function of the anisotropy parameter. 
Figure 2: The ground state energy as the function of the anisotropy parameter. 
Figure 3: The spontaneous magnetization as the function of the anisotropy parameter. 
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